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SYMBOLS 
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CCL  -  Convective  condensation  level 

Cp  -  Heat  capacity  of  air  at  constant  pressure 

CT  -  Convective  temperature 

E  -  Actual  vapor  pressure 

ES  -  Saturation  vapor  pressure 

FR  -  Relative  humidity 

I ,  J  ,k  -  I  ndexes 

L  -  Latent  heat  of  vaporization  of  water 
LCL  -  Lifting  condensation  level 
M  -  Saturation  vapor  pressure  over  water 
P  -  Pressure 

P*  -  Pressure  correction 

PB  -  Pressure  at  the  bottom  of  a  layer 

PC  -  Pressure  at  the  convective  condensation  level 

PI  -  Pressure  at  the  intersection 

PM  -  Pressure  at  the  top  of  the  mixing  layer 

PS  -  Surface  pressure 

PT  -  Pressure  at  the  top  of  a  layer 

R  -  Gas  constant 

T  -  Temperature 

T*  -  Temperature  correction 

TD  -  Dewpoint  temperature 

TDS  -  Dewpoint  temperature  at  the  surface 

Preceding  page  blank 
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T  ^  -  Temperature  on  a  dry  adiabatic  curve 
TE  -  Psuedo  equivalent  temperature 
Tl  -  Temperature  at  an  Intersection 

T^  -  Temperature  at  the  top  of  the  mixing  layer 

T  -  Temperature  on  a  mixing  ratio  curve 
MR 

T^  •  Temperature  on  a  saturation  adlabat  curve 

TW  -  Wet  bulb  temperature 

W  -  Mixing  ratio 

W  -  Mean  mixing  ratio 

X  -  Coordinate 

Y  -  Coord i nate 

Z  -  Thickness  of  a  layer 

0  -  Potential  temperature 

0E  -  Equivalent  potential  temperature 

0^  -  Parameter  for  saturation  adlabat 
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INTRODUCTION 


The  increasing  availability  of  computing  facilities,  programmable 
calculators,  and  minicomputers  allows  many  of  the  computations 
currently  performed  by  manual  graphics  to  be  done  by  computer. 

This  paper  discusses  numerical  methods  of  computing  meteorological 
quantities  which  are  usually  manually  derived  from  analysis  on  a 
SKEW-T,  log  p  DIAGRAM  (or  SKEW-T).  The  numerical  methods  used  were 
selected  for  their  simplicity  and  accuracy.  A  mathematical  char¬ 
acterization  of  the  SKEW-T  and  algorithms  for  computing  several 
meteorological  quantities  are  presented.  Finally,  a  discussion  of 
the  application  of  these  methods  and  a  FORTRAN  program  listing  to 
accomplish  the  computations  are  included. 

THE  SKEW-T,  log  p  DIAGRAM 

The  SKEW-T,  log  p  DIAGRAM  CG  is  a  thermodynamic  chart  with  five 
families  of  curves,  five  types  of  scales,  and  three  data  blocks. 

Various  numerical  Information  is  also  tabulated  on  the  SKEW-T.  This 
paper  is  concerned  chiefly  with  the  five  families  of  curves  which  are 
associated  with  pressure,  temperature,  dry  adiabat,  saturation 
adiabat,  and  mixing  ratio. 

The  first  two  families  of  curves,  temperature  and  pressure,  are  used 
to  locate  points  on  the  chart.  An  arbitrary  coordinate  system  has  been 
selected  to  measure  distances.  Let  the  origin  correspond  to  the  point 
at  a  temperature  of  OC  (centigrade)  and  a  pressure  of  1000  mb  (mil¬ 
libars).  Take  the  X  direction  to  be  parallel  to  the  pressure  lines 
(horizontal),  with  positive  X  to  the  right.  The  point  at  a  temperature 
of  1C  and  a  pressure  of  1000  mb  is  on  the  positive  side  of  the  origin. 
The  Y  direction  is  perpendicular  to  the  X  direction.  Positive  Y  is 
towards  lower  pressures  (up).  A  point  on  the  chart  which  is  specified 
by  its  temperature  and  pressure  may  be  transformed  to  X,Y  coordinates 
by  Eqs.  (I)  and  (2).  The  components  in  the  X,Y  coordinate  system  are 
given  in  inches. 

X  =  .  I  408T  -  10.53975  log|Q  P  +  31.61923  (!) 

Y  =  -11.5  log ! 0  P  +  34.5  (2) 


The  X,Y  coordinates  have  been  scaled  to  USAF  SKEW-T,  log  p  DIAGRAM 
DOD-WPC-9- 16- I .  See  [l]. 


The  remaining  three  families  of  curves,  TW,  TSA,  and  TDA,  are 
given  In  Table  I.  The  temperatures  are  specified  as  a  function  of 
pressure  and  a  parameter,  the  parameter  serving  as  a  means  of 
specifying  a  particular  curve  of  the  family. 

The  temperature  T  at  an  arbitrary  pressure  on  a  saturation  adiabat 
Is  determined  by  the  bisection  method. r  The  temperature  is  assumed 
to  lie  In  the  range  -80C  to  40C.  An  Initial  guess  of  -20C  is 
made  anj  the  correction,  T*,  computed.  The  correction  term  decreases 
by  a  fzctor  of  1/2  after  each  correction.  Terminating  after  13 
corrections  gave  satisfactory  results.  The  algorithm  for  computing 
the  temperature  on  a  saturation  adiabat  Is  based  on  Eq.  (3): 


0 


( 6E )  •  EXP(- 


(3) 


The  latent  heat  of  vaporization  L  and  the  heat  capacity  of  air 
at  constant  pressure  Cp,  are  considered  constant.  Since  It 
is  not  known  how  accurately  the  saturation  adiabat  could  be 
determined  from  Eq .  (3),  Table  2  was  constructed  using  List 
C23  as  a  standard.  The  temperature  on  an  arbitrary  mixing  ratio 
curve  W  is  computed  by  first  computing  the  saturation  vapor 
pressure  M.  An  approximation  to  the  inverse  saturation  vapor 
pressure  function  is  then  used  to  compute  the  temperature. 

In  addition  to  the  algorithms  which  generate  the  curves  for 
each  family,  it  is  necessary  to  have  algorithms  which  determine 
which  curve  in  a  family  passes  through  an  arbitrary  point  (T,P). 
Algorithms  to  accompl ish  this  are  given  in  Table  3. 


■J-The  bisection  method  Is  a  numerical  technique  which  decreases  the 
difference  between  the  upper  and  lower  estimates  by  a  factor  of 
1/2  per  iteration. 
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TABLE  2 


TEMPERATURE  AND  ERROR  ON  SELECTED  SATURATION 
ADIABATS  AT  SELECTED  PRESSURES 


Pressure 

(mb) 

Temperature 

(C) 

Error 

(C) 

1000.0 

40.0000 

70!  .5 

29.9877 

.0122 

490.7 

19.9536 

.0463 

344.7 

9.9194 

.0805 

245.4 

-.1733 

-.1733 

179.6 

-10.2221 

-.2221 

1000.0 

30.0000 

733.0 

19.9829 

.0170 

544.0 

9.9633 

.0366 

412.4 

-.0561 

-.0561 

321.4 

-10.0756 

-.0756 

257.7 

-20.1538 

-.1538 

212.0 

-30.2612 

-.2612 

177.6 

-40.3247 

-.3247 

1000.0 

20.0000 

770.0 

9.9780 

.0219 

606.0 

-.0561  ■ 

-.0561 

489.0 

-10.0463 

-.0463 

403.0 

-20.1245 

-.1245 

338.0 

-30.1879 

-.1879 

286.4 

-40.2368 

-.2368 

243.5 

-50.2709 

-.2709 

206.8 

-60.2612 

-.2612 

174.7 

-70.2661 

-.2661 

1000.0 

10.0000 

805.0 

-.0415 

-.0415 

663.0 

-9.9877 

.0122 

554.0 

-20.0952 

-.0952 

470.0 

-30.0268 

-.0268 

400.0 

-40.1196 

-.1196 

341.0 

-50.1538 

-.1538 

289.9 

-60.1586 

-.1586 

245.1 

-70.1489 

-.1489 

205.7 

-80.1098 

-.1098 

171.0 

-90.1 147 

-.1147 
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TABLE  2  (con.) 


Pressure 

(mb) 

Temperature 

(C) 

Error 

(C) 

1000.0 

.0000 

833.0 

-9.9731 

.0268 

703.0 

-19.9926 

.0073 

599.0 

-29.9829 

.0170 

51  1  .0 

-40.1 196 

-.1196 

436.4 

-50.1391 

-.1391 

371.3 

-60.1293 

-.1293 

314.0 

-70.1 196 

- . 1 1 9b 

263.5 

-80.0952 

-.0952 

219.1 

-90.0854 

-.0854 

1000.0 

-10.0000 

849.0 

-20.0073 

-.0073 

726.0 

-29.9829 

.0170 

621  .0 

-40.0756 

-.0756 

531.2 

-50.0512 

-.0512 

452.2 

-60.0415 

-.0415 

382.4 

-70.0463 

-.0463 

266.9 

-89.9975 

-.0024 

1000.0 

-20.0000 

856.8 

-30.0122 

-.0122 

734.8 

-40.0170 

-.0170 

628.0 

-50.0366 

-.0366 

535.3 

-60.0268 

-.0268 

452.8 

-70.0170 

-.0170 

380.0 

-80.0073 

-.0073 

316.0 

-89.9829 

.0170 
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TABLE  3 


DETERMINING  A  CURVE  THROUGH  A  GIVEN  POINT 


FAMILY 

PARANETER  "FOR  CURVE  PASSING  “THROUGH  (f,P) 

Dry  Adiabat 

e  =  T{'000).28b 

MI x I ng  Ratio 

...  622  ESAT(T) 

P  -  ESAT (T) 

Saturation  Adiabat 

T. 1000  ..288 
v  p 

EXP  (  bWjT'P)) 

b  =  -2.6518986 

NOTE:  T  Is  In  Kelvin.  K  *  C  +  273.16  (see  Table  I  for  a  definition 
of  ESAT) 
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ALGORITHMS  FOR  SELECTED  METEOROLOGICAL  QUANTITIES 


Several  meteorological  quantities  which  are  usually  manually 
derived  from  an  analysis  of  a  SKEW-T  were  selected  for  discussion. 
Algorithms  are  presented  for  corrputing  these  meteorological 
quantities.  The  selection  of  symbols  Is  somewhat  different 
than  is  customary  because  of  current  symbol  limitations  on 
corrputers.  But  by  referring  to  the  List  of  Symbols,  the  reader 
will  have  no  difficulty.  Units  are  the  same  as  those  used 
on  the  SKEW-T  [l]. 


Mi xl ng  Ratio:  W 

The  mixing  ratio  W  Is  computed  from  the  pressure  P  and  the 
dewpoint  temperature  TD  by  using  the  function  ESAT,  which  is 
def i ned  in  Table  I . 


.  622  ESAT(TD) 

w  "  nw TTDT 


(4) 


TD  i  s  in  degrees  Ke  I  vi  n,  the  pressure  P  in  mi  1 1  i  bars ,  and  W  in 
grams  of  water  vapor  per  kilogram  of  dry  air.  The  saturat  ion 
mix!  ng  rat  io  is  obtained  by  using  the  dry  bulb  temperature 
in  place  of  the  dewpoint  temperature. 

Relative  Humidity:  FR 

The  relative  humidity  is  computed  from  the  temperature  T  and 
the  dewpoint  temperature  TD  by  using  ESAT.  Both  T  and  TD 
are  in  degrees  Kelvin. 

FR  =  100  (ESAT(TD)/ ESAT(T) )  (5) 

Saturation  Vapor  Pressure:  ES 

ESAT  gives  the  saturation  vapor  pressure  in  millibars  from  the 
dry  bulb  temperature  T,  which  is  in  degrees  Kelvin. 


ES  =  ESAT(T) 


(6) 
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Actual  Vapor  Pressure:  E 

The  dewpoint  temperature  TD  Is  used  In  place  of  T  In  Eq.  (6). 
Potential  Temperature:  0 

The  potential  temperature  Is  computed  from  the  dry  bulb  temperature 
T  in  Kelvin  and  the  pressure  P  In  millibars. 


e  .  T  -288 


(7) 


The  Wet  Bulb  Temperature  and  Wet 
Bulb  Potential  Temperature:  TW,0W 

The  wet  bulb  temperature  Is  approximated  by  calculating  the  psuedo 
wet  bulb  temperature.  The  arguments  are  surface  dewpoint  temper¬ 
ature,  surface  temperature,  and  pressure,  which  are  symbolized 
by  TDS,  TS,  and  PS,  respectively.  TOS  and  TS  are  In  Kelvin, 

P  in  millibars.  First  a  mixing  ratio  curve  W,  which  passes 
through  (TDS,  PS),  Is  determined.  From  Table  3  we  have: 


622  ESAT  (TDS ) 
P  -  ESAT  (TDS) 


(8) 


Next  a  dry  adiabat,  which  passes  through  (TS,  PS),  Is  determined. 
Again  by  referring  to  Table  3  we  have: 


0  =  TS(- 


PS  . .288 


1000 


(9) 


Two  curves  have  now  been  specified:  T^  (W,  P)  and  T».  (0,  P). 
The  next  step  Is  to  locate  the  pressure  at  which  the  curves 
Intersect.  This  is  done  by  an  iterative  procedure.  An  initial 
guess  that  the  Intersection  pressure  PI  Is  equal  to  the  surface 
pressure  is  made.  A  correction  Is  computed  and  a  revised  guess 
made.  When  (TUCJ  -  T  )2  Is  loss  than  .0001  degrees,  the  process 
Is  terminated. 


PI,  =  PS 


(10) 
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+  p* 


(11) 


PI,  =  PI 


2(-02(TNK(W'FV 


TDA(6'Pk))) 


(12) 


It  is  found  that  six  iterations  were  sufficient  to  compute  PI 
to  within  I  mb.  Once  the  pressure  and  hence  temperature  at 
the  Intersection  are  known,  a  saturation  adiabat  through  the 
intersection  point  (Tl ,  PI)  is  found.  Referring  to  Table  3  we 
have: 


, 1 000,. 288 
=  TI  ~PT~~ 

exp  (Mllj-fil, 


(13) 


Finally  by  fol  low i ng  th i s  saturation  adiabat  to  the  surface 
pressure  PS  and  to  1000  mb,  we  get  the  wet  bulb  temperature 
TW  and  the  wet  bulb  potential  temperature  GW,  respectively. 


Tw  =  tsa  (eS' 

PS) 

(14) 

0W  =  TSA  (0S' 

1000) 

(IS) 

The  Psuedo  Wet  Bulb  Temperature 
and  Pseudo  Wet  Bulb  Potential  Temperature: 

TPW,  6PW 

Refer  to  the  wet  bulb  temperature  and  wet  bulb  potential  temperature 
above . 


The  Equivalent  Potential  Temperature:  0E 

The  equivalent  potential  temperature  is  computed  from  the  same 
quantities  used  to  compute  the  wet  bulb  temperature,  i.e.,  the  surface 
pressure,  dewpoint  temperature,  and  actual  temperature.  First 
compute  the  wet  bulb  temperature  TW.  The.  equ i va lent  potential 
temperature  can  then  be  computed  by  the  same  process  used  to 
determine  the  parameter  9_  of  a  saturation  adiabat  through  (TW,  PS). 
Referring  to  Table  3  we  h3ve: 


14 


6E  = 


(16) 


.1000.. 288 
TW  PS 

EXP 


The  Psuedo  Equivalent  Temperature:  TE 

First  the  equivalent  potential  temperature  6E  Is  computed.  The 
psuedo  equivalent  temperature  is  then  given  by 

Tc  __  ,  PS  ,.288 

TE  '  eE  lTooo  1  U‘ 

Thickness  of  a  Layer:  Z 

It  is  assumed  that  the  temperature  and  dewpoint  temperature  are  known 
at  N  distinct,  decreasing  pressures.  Thicknesses  are  computed 
In  meters  from  the  surface.  The  trapezoidal  rule  Is  used  to 
I  ntegrate 


i  nr  i 

Z  =  y|g  J  [T  +  . 6078*W*T/( 1 000  +  W)]d  InP 


See  Table  I  for  a  definition  of  W(T,  P). 

Rewriting  Eq .  (18)  and  noting  that  W<<I000  gives  Eq.  (19),  which 
is  used  to  perform  the  computation  of  Z: 

(T  +  T  +  6.078-  10“6*  (W  -T  +  W  -T  )) 

Z  =  29.2857  — - - - - - - - - - — - -  •  ln(P  /P  ) 

2 


(1  +  T  +  6.078-10  -(W,-T  +  W  •  T  )) 

+  - - - - - - - - - - - ln(P?/P  ) 

2  ^ 


+  . . .  + 


(T  +  T  .  +  6.078-10  U-(W  -T  +  W  ..-T  ..)) 
n  n+ I  n  n  n+l  n+ I 


•T  +  W  .  - T  A. )) 

^ — - a*1-"*1  .  ln(P  /P  ,  )  (13) 

n  n+ 1 


The  Lifting  Condensation  Level:  LCL 

The  lifting  condensation  level  is  computed  In  the  same  manner 
that  PI  was  computed  for  the  wet  bulb  temperature.  Eqs.  (8),  (9), 
(10),  (II),  and  (12)  are  used.  (Tl,  PI)  locate  the  LCL. 

The  Convective  Condensation  Level:  CCL 

It  Is  assumed  that  the  temperature  and  dewpoint  temperature  are 
known  at  N  distinct,  decreasing  pressures.  The  pressure  at  the 
top  of  the  mixing  layer  PM  must  be  greater  than  P  ,  the  last 
pressure  level.  Since  PM  is  bounded  by  P.  and  P  ,  there  Is 
a  K  such  that  n 


Pk  >  PM  >  PRf 


(20) 


First  the  mean  mixing  ratio,  W,  in  the  P  -  PM  layer  is  computed: 

k- 

w.  lJ| 


*  [W‘V  pi’  +  “(Tm-  fW]  <lnP;  -  lnPm> 


2(  In  P  -  In  PM) 


W(T.  ,  PJ  +  W(T  »  pM)  (In  P,  -  In  PM) 
k  k _ m _ _ k 

2( In  Pj  -  In  PM) 


(23) 


The  intersection  of  T(^  (W,  P)  and  the  curve  defined  by 


Ts  <P)  ■  Tk  - 


(Vi  ~  Tkliinp  -  iniy 

<lnPk  -  lnpk+l> 


(22) 
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(k  is  chosen  such  that  P  >  P  _>  P.  . )  defines  the  convective 
condensation  level.  This  i ntersecTlon  can  be  found  first 
systematically  comparing  the  difference  between  (W,  P.) 
and  Tj.(P. )  until  the  smallest  i  is  found,  such  that 


Tf*(W,  P,)  -  Ts  <Pj) 


<  0 


and 


T,* 


pm> 


-  T. 


iP|+l> 


(23) 


(24) 


A  bisection  msthod  Is  used  to  determine  PC,  the  pressure  at  the 
CCL.  An  initial  guess  PC.  is  made,  tested  to  see  If  T^„(W,  PC  ) 
equals  ( PC  ^ ) ,  and  if  nor  ,  corrected. 


PC,  =  .5  (P.  + 


pm> 


(25) 


PC.  =  PC.  ,  +  P*  (corrector) 
J  J  '  J“  • 


(26) 


pk  -  ~h—  SIGN  <TW)  (w'  PV  ■  Ts  (Pk,) 


(27) 


Ten  corrections  are  made. 

The  Convective  Temperature:  CT 

First  the  pressure  PC  at  the  convective  condensation  level  is 
computed.  The  temperature  at  the  CCL,  TC,  is  computed  from  PC 
and  W: 


TC  =  Tf* 


(W,  PC) 


(28) 
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A  dry  adiabat  is  determined. 


6 


1000. .288 

PC 


(29) 


Finally,  the  convective  temperature  CT  is  corrputed  from  6  and  the 
surface  pressure  PS: 


CT  =  Q  ( 


PS 

1000 


j.288 


(30) 


APPLICATIONS 

The  algorithms  are  useful  for  data  reduction  purposes.  The  memory 
and  speed  requirements  are  not  excessive  and  most  computations 
can  be  carried  out  satisfactorily  on  a  programmable  calculator. 

In  addition  to  data  analysis,  the  algorithms  are  useful  for 
generating  backgrounds  for  the  presentation  of  data.  An  example 
of  a  computer  generated  background  and  plotted  sounding  is  given 
in  Fig.  I.  Computation  of  selected  meteorological  quantities  from 
the  sounding  in  Fig.  I  are  presented  in  Table  4.  A  table  of  CCL 
temperatures,  pressures,  and  heights  was  computed  using  an 
arbitrary  decrement  of  -25  mb  for  the  pressure  at  the  top  of  the 
mixing  ratio. 
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S  I f  IS? 


100 
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300 

400 


SKEW  T.LOC  P 


Figure  I.  A  computer-generated  SKEW-T  background  with  a  hypothetical 

profile  of  temperature  and  dewpoint  temperature.  Horizontal 
I  int.*  are  pressure.  Positive  sloped  lines  are  temperature 
(solia)  and  mixing  ratio  (dashed).  Negative  sloped  curves 
are  dry  adiabat  (dashed)  and  saturated  adiabat  (solid). 
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TABLE  4 


AN  EXAMPLE  OF  A  VERTICAL  SOUNDING 


Press. 

1013 

953  950 

942 

920 

843  777 

745 

691 

620  333 

210 

Height 

0 

536  554 

626 

827 

1553  2213 

2549 

3145 

3984  8664 

12047 

Pot.  T. 

19.3 

22.3  23.3 

23.2 

21  .3 

19.8  19.8 

24.0 

24.5 

26.2  74.1 

115.0 

Temp . 

20.4 

18.2  19.0 

18.2 

14.4 

5.8  -0.7 

-0. 1 

-5.5 

-12.3  -20.1 

-25.5 

Dew  pt. 

18.2 

14.4  6.0 

-0.8 

-0.6 

-5.2  -12.7 

-20.1 

-25.5 

-30.3  -28.1 

-32.5 

R.H. 

87 

79  43 

28 

36 

45  40 

21 

19 

21  49 

52 

Mi x  Rat lo 

13. 1 

10.9  6.09 

3.86 

4.10 

3.01  1.91 

1  .04 

.70 

.51  1.14 

1.17 

Sat.  V.P. 

24.0 

20.9  22.0 

20.9 

16.4 

9.2  5.8 

6.0 

4.1 

2.4  1.2 

0.8 

V.  Press. 

20.9 

16.4  9.2 

5.8 

6.0 

4.1  2.4 

1  .2 

0.8 

0.5  0.6 

0.4 

W.  Bulb 

Temp 

19.0 

15.9  11.7 

8.8 

7.1 

1.2  -4.7 

-5.9  ■ 

-10.2 

-15.7  -22.1 

-27.4 

Pot.  W.B. 

Temp 

Equlv.  Pot. 

18.5 

17.8  13.7 

1 1  .4 

10.9 

9.0  7.7 

8.4 

8.2 

8.7  23.9 

31  .1 

Temp 

56.7 

54.2  41.3 

34.8 

33.5 

28.9  25.8 

27.3 

26.9 

28.0  78.6 

120.3 

LCL  at 

Temp . 

17.9  Press. 

983 

Height  260 

Mixing  Layer 

CCL 

Convect 1 ve 

Temp 

Mean  Mixing 

Ratio 

Press 

Height 

Press 

Height 

988 

215 

931 

725 

23.4 

12.64 

963 

437 

928 

757 

23.  1 

12.19 

938 

662 

91 1 

908 

22.3 

10.68 

913 

890 

882 

1173 

21  .8 

8.95 

888 

1121 

863 

1361 

21  .4 

7.86 

863 

1359 

847 

151  I 

21  .1 

7.08 

838 

1602 

832 

1660 

21  .0 

6.47 

813 

1846 

817 

1807 

21  .0 

5.95 

788 

2099 

803 

JQ/7 

21  .0 

5.50 

763 

2358 

790 

2075 

21  .0 

5.08 

Units:  temperature  C;  pressure  millibar;  mixing  ratio  grams/kl  logram;  height  meters. 
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